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Abstract New exact analytic solutions are introduced for the rotational 
motion of a rigid body having two equal principal moments of inertia and 
subjected to an external torque which is constant in magnitude. In partic- 
ular, the solutions are obtained for the following cases: (1) Torque parallel 
to the symmetry axis and arbitrary initial angular velocity; (2) Torque per- 
pendicular to the symmetry axis and such that the torque is rotating at a 
constant rate about the symmetry axis, and arbitrary initial angular veloc- 
ity; (3) Torque and initial angular velocity perpendicular to the symmetry 
axis, with the torque being fixed with the body. In addition to the solutions 
for these three forced cases, an original solution is introduced for the case 
of torque-free motion, which is simpler than the classical solution as regards 
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its derivation and uses the rotation matrix in order to describe the body ori- 
entation. This paper builds upon the recently discovered exact solution for 
the motion of a rigid body with a spherical ellipsoid of inertia. In particular, 
by following Hestenes' theory, the rotational motion of an axially symmetric 
rigid body is seen at any instant in time as the combination of the motion 
of a "virtual" spherical body with respect to the inertial frame and the mo- 
tion of the axially symmetric body with respect to this "virtual" body. The 
kinematic solutions are presented in terms of the rotation matrix. The newly 
found exact analytic solutions are valid for any motion time length and rota- 
tion amplitude. The present paper adds further elements to the small set of 
special cases for which an exact solution of the rotational motion of a rigid 
body exists. 

Keywords rigid body dynamics, kinematics, rotation, integrable cases of 
motion, revolution ellipsoid of inertia, analytical reduction 

1 Introduction 

The study of angular motion of a rigid body is at the very basis of clas- 
sical physics. Besides the theoretical interest, the subject is of considerable 
practical importance for the fields of astronautics and celestial mechanics. 

The problem of the rotational motion of a rigid body can be divided into 
two parts. The dynamic problem aims to obtain the angular velocity of the 
body with respect to an inertial reference by starting from the knowledge of 
the initial angular velocity and the history of the applied torque. On the other 
hand, the kinematic problem focuses on determining the current orientation 
of the body from the knowledge of the initial orientation and the history of 
the angular velocity. 

The exact analytic solution for the rotational motion of a rigid body exists 
only for a small number of special cases. 
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When no external torque (other than the one due to gravity) acts on the 
body, exact solutions exist for both the dynamic and kinematic problems 
only in the following three cases: 

1. Eulcr-Poinsot case: torque- free rotation of an asymmetric rigid body with 
a fixed point coincident with the center of mass (Euler 1758; Poinsot 1834; 
Jacobi 1849; Leimanis 1965; Morton et al. 1974); 

2. Lagrange-Poisson heavy top case: rigid body under gravity force with 
two equal principal moments of inertia at the fixed point and the center 
of mass along the third axis of inertia (Lagrange 1788; Poisson 1813; 
Golubev 1960; Leimanis 1965); 

3. Kovalevskaya heavy top case: rigid body under gravity force with two 
equal principal moments of inertia at the fixed point, value of the third 
moment of inertia equal to half of the value of the moment inertia about 
the other two axes, and the center of mass in the plane of equal moments 
of inertia (Kovalevskaya 1890; Leimanis 1965). 

On the contrary, when the rigid body is subjected to an external torque, 
the complete (i.e. for both the kinematics and the dynamics) exact analytic 
solution exists only for a rigid body with spherical ellipsoid of inertia sub- 
jected to a constant torque and arbitrary initial angular velocity (Romano 
2008), in addition to the straightforward single degree-of- freedom case of 
rotation with a torque along a principal axis of inertia and initial angular 
velocity along the same axis. 

Furthermore, a partial (limited to the dynamic problem) exact analytic 
solution exists for the case of an axially symmetric body subjected to a 
constant torque (Longuski 1991, Tsiotras and Longuski 1991). In particular, 
Longuski 1991 and Tsiotras and Longuski 1991 use a complex form expression 
of the Euler's dynamic equations and give an exact solution for the angular 
velocity involving the Fresncl integral function, and an approximate solution 
for the kinematic problem in terms of Euler angles. Nevertheless, no exact 
analytic solution has been found for the kinematic problem. 
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Several researchers have proposed approximate solutions. In particular, 
Longuski 1991, and Longuski and Tsiotras 1995 use asymptotic and series 
expansion to get the approximate solution of the dynamics and kinematics 
for both the constant torque case and the more difficult case of time- varying 
torque. Livneh and Wie 1997 analyze qualitatively, in the phase space, the 
dynamics of a rigid-body subjected to constant torque. As regards the kine- 
matic problem, Iserles and N0rsett 1999 study the solution in terms of series 
expansion for the more general problem of solving linear differential equation 
in Lie Groups, building upon the work of Magnus 1954. Finally, Celledoni 
and Saefstroem 2006 propose ad-hoc numerical integration algorithms. 

For the first time, to the knowledge of the author, the present paper 
introduces the complete exact analytic solutions for the rotation of an axially 
symmetric rigid body in the following three cases 

1. External torque constant in magnitude and parallel to the symmetry axis, 
with arbitrary initial angular velocity and body orientation. 

2. External torque constant in magnitude, perpendicular to the symmetry 
axis and rotating at a constant rate about the symmetry axis, with arbi- 
trary initial angular velocity and body orientation. 

3. External torque and initial angular velocity perpendicular to the symme- 
try axis, with the torque constant in magnitude and fixed with the body 
and arbitrary initial body orientation. 

The solution for the kinematics are given in terms of the rotation matrix. 

This paper builds upon the recently found exact analytic solution for 
the motion of a rigid body with a spherical ellipsoid of inertia (Romano 
2008), and the Hestenes' method of decomposing the rotational motion of 
an axially symmetric rigid body into the combination of the motion of a 
"virtual" spherical body with respect to the inertial frame and the motion 
of the axially symmetric body with respect to this "virtual" body (adapted 
from Hestenes 1999). 
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In addition to the three cases listed above, an original complete solution 
of the Eulcr-Poinsot's classical case of torque- free motion is also presented 
in this paper, which is simpler than the classical solution (Leimanis 1965) as 
regards its derivation and uses the rotation matrix in order to describe the 
body orientation. 

The newly introduced exact analytic solutions are of high theoretical in- 
terest, as they add new elements to the set of few special cases for which 
an exact solution of the rotational motion of a rigid body exists. Addition- 
ally, they are also of interest from the applied mathematics point of view as 
they constitute significant new comparison cases for the validation and error 
analysis of approximate algorithms. Indeed, the analytic solutions presented 
in this paper arc valid for any length of time and rotation amplitude: they 
only require the numerical evaluation of mathematical expressions, without 
needing any numerical propagation. 

Finally, the solutions introduced in this paper are valid for any rigid body 
which has an ellipsoid of inertia of revolution, i.e. two equal principal mo- 
ments of inertia. This class of bodies, having a dynamic axial symmetry, 
contains as a subset the class of bodies having geometric axial symmetry 
(Hestenes 1999). However, following the common practice adopted in liter- 
ature, in this paper the definition of "axially symmetric body" and "body 
with revolution ellipsoid of inertia" are used interchangeably. 

The paper is organized as follows: Section 2 briefly introduces the dynamic 
and kinematic equations for the rotational motion of a generic rigid body and 
Section 3 summarizes the known results for a rigid body having spherical 
ellipsoid of inertia. Section 4 introduces the method used to analyze the 
motion of a rigid body with revolution ellipsoid of inertia, and Section 5 
presents the new exact analytic solutions listed above. Finally, Section 6 
concludes the paper. 
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2 General Rigid-Body Problem Statement 

For a generic rigid body the Eulcr's rotational equation of motion, in vectorial 
form, is (Goldstein 1980) 

h = m (1) 

where the the dot symbol denotes the time derivative with respect to an 
inertial reference frame, m is the resultant external torque acting on the 
body, and 

h=Icj (2) 

is the absolute angular momentum of the body, I is the inertia dyadic of the 
body with respect to its center of mass, and lo is the angular velocity of the 
body with respect to an inertial reference frame. 

By resolving all of the vectors and the inertia dyadic along a body-fixed 
Cartesian coordinate system B with axis equal to the principal axes of inertia, 
Eq. 1 can be written in scalar form as 

hp = (h -h)qr + mi 

hq = (h - h) r p + m 2 (3) 
hf = (h - h)pq + m 3 , 

where {Ii : i = 1, 2, 3} are the principal moments of inertia, {p, q, r} are the 
components in the coordinate system B of u, and elements of the column 
matrix b ojbn- Finally {m t : i = 1,2,3} are the components of m along the 
same axes. 

The rotation matrix Rnb G SO (3) from the body fixed coordinate sys- 
tem B to an inertial coordinate system N obeys the following differential 
equation (Leimanis 1965) 



(4) 
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where 
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In general, Equations 3 and 4, which solve the dynamic and kinematic motion 
problem of a rigid body respectively, do not have exact analytic solutions. In 
particular, the following straightforward direct solution of Eq. 4 



where exp(») indicates the matrix exponential, is only valid for the very par- 
ticular cases when the matrix Q(t) commutes with its integral (Adrianova 
1995). This commutativity is lost when the initial angular velocity is arbi- 
trary; therefore, Eq. 6 is not valid anymore in that case. The limitation in 
the validity of the solution of Eq. 6 is not clearly stated in Bodewadt 1952 
and Leimanis 1965, as already noticed in Longuski 1984. 

3 Exact analytic solution for the motion of a rigid body with a 
spherical ellipsoid of inertia 

This section summarizes the results introduced by Romano 2008 for the 
motion of a rigid body with spherical ellipsoid of inertia and subjected to 
a constant torque. The developments reported in this section will be used 
in the following sections of this paper in order to analyze the motion of an 
axially symmetric rigid body. 

For a rigid body having a spherical ellipsoid of inertia, the kinematic 
differential equations in terms of the stereographic complex rotation variables 
Wk (Tsiotras and Longuski 1991, Schaub and Junkins 2003) are (Romano 
2008) 





1 

Wk = - 
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Theorem 1 Given pQ,q a ,r and U real numbers, with p and q not both 
zero, the general solution for each one of the Eq. 7, governing the rotational 
kinematics of a rigid body with spherical ellipsoid of inertia, initial angular 
velocity components p(0) =Po, q(0) = qo, and r(0) = ro along the three body 
fixed axes o~\,Gi and a^, and subjected to a constant torque U, normalized 
by the value of the moment of inertia and directed along the axis 03, is the 
following 



with 



w(t,c)= [\ + i)VU , [6z + G(z,c)}, (8) 



(9) 



, 2 1F1 (g^, f ; z 2 ) (u - 1) z 2 + 6 1F1 (1 - § , |; z 2 ) cvz - 3 1 F 1 

where 1F1 denotes the confluent hypergeometric function (Lebedev 1965), c E 
C is £/ie constant of integration and 

z . = (l + i)(r + Ut) 1 »(Po+go) flQl 

2^ ' 417 ' 1 J 



Corollary 1 The solution in terms of the rotation matrix, which corresponds 
to the solution given by theorem 1 in terms of stereographic rotation variables, 
is 

R N B(t) = [r kj (t)], k, 3 = 1,2,3, (11) 

with 

i(w k -Wk) w k +Wk~ l-\w k \ 2 , 10 _ 

l + |w fe | 2 l + \w k \ A l + \w k \ 2 

(12) 

where w k — w(t,c k ), being w(t,c k ) given by Eq. 8 with the following values 
for the initial conditions (obtained by considering the rotation matrix Rnb 
to be equal to the identity matrix at the initial time t = 0, without loss of 
generality). 

(1 + i) f i p i (^>bw) [6r 2 + 3(po-iq )r a + 6iU]+2g 



Cl 



6\/C7 | 2 !Fx (l - §, §; g) ^r + iF a |; g) (p - z<z + 2r ) 
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2 ' 2 ' 2(7 

where v is defined as in Eq. 10, and 

'3 








^ 5 irg 



(14) 



2 '2' 217, 

In the general case when the constant torque is not directed along the axis 
(T3 , and the coordinate system B is not coincident with N at the initial time, 
the results of Theorem 1 and Corollary 1 are still applicable, as reported in 
the following corollary, which constitutes a new development with respect to 
Romano 2008. 

Corollary 2 Let us consider that the direction of the external torque is iden- 
tified in the coordinate system B by the unit vector 



> - 



(15) 



and that the initial orientation of B with respect to the inertial coordinate 
system N is given by Rbn{0)- Then, by exploiting the property of successive 
rotations, it results 



RbnW — Rbk RKN{t) Rkb Rbn(0), 



(16) 



where Rkb is any time-independent rotation matrix which brings the direc- 
tion of the torque to coincide with the third axis of the auxiliary body-fixed 
coordinate system K . In other words, the matrix Rkb needs to satisfy the 
condition 



(17) 





Ul 






Rkb < 
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For instance, Rkb can be obtained by combining a first elementary rotation 
of an angle a = arctan 2(i*2, Mi) (indicating by arctan2 the four quadrant 
inverse tangent) about the third axis of B with a second elementary rotation of 
an angle (3 = arctan 2(u2/sin(a), U3) about the resulting second axis, yielding 



cos (a) sin (a) 
— sin (a) cos (a) 
1 



(18) 



cos (/3) - sin ((3) 
Rkb = 10 

sin (j3) cos (f3) 

Furthermore, Ri<N(t) is obtained by transposing the resulting matrix of Eq. 11 
of Corollary 1, with the subindex B substituted by K , and the following val- 
ues of the components along the coordinate system K of the initial angular 
velocity of the body with respect to the inertial frame 



(19) 



where the prime symbol is used to denote the initial conditions expressed in 
the B coordinate system. 



H 








> — Rkb ' 











4 Analysis of the rotational motion of a rigid body with a 
revolution ellipsoid of inertia 

Let us assume, from now on, that the rigid body has ellipsoid of inertia 
which is a revolution ellipsoid. Consequently, let us assume that the principal 
moments of inertia satisfy the following relation 

h=h=I± h- (20) 

By following the development of Hestenes 1999, the absolute angular mo- 
mentum of the body can be expressed by 

h=Xu = Iw + (h-I)(u-e)e, (21) 

where e is the unit vector along the direction of the symmetry axis. 
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Additionally, the angular velocity can be seen, at any instant of time, as 
the sum of two vectorial components, whose one is parallel to the angular 
momentum vector (tv h ) and the other one is parallel to the body symmetry 
axis (o; e ), i.e. 

w = w fe +w e . (22) 
In particular, from Eq. 21, it yields 

h (I-I 3 ) fh 

u h = j, u e = 



j^l l=-e) e = A(w h -e)e, (23) 



with the constant A defined as A = (I — I^/I^. 
From Eq. 23 it results 

h = Iu, h , (24) 
and, from the Euler's equation (Eq. 1), 

h = Iuj_ h = m (25) 

Equations 22 and 25 are the mathematical expression of the Reduction 
Theorem (Hestenes 1999) which can be formally stated as follows: the evo- 
lution in time of the absolute angular momentum of an axially symmetric 
body subjected to an external torque vector m is corresponding to the mo- 
tion of a "virtual" homogeneous spherical body with the same value of the 
transversal inertia of the axially symmetric body and subjected to the same 
external torque. 

In other words, to exemplify, the rotational motion of an axially symmet- 
ric rigid body (with two equal moments of inertia / and third moment of 
inertia / 3 ) is decomposed into the combination of the motion of a "virtual" 
spherical body (with principal moment of inertia /) with respect to the iner- 
tial frame and the spinning motion, about its symmetry axis, of the axially 
symmetric rigid body with respect to this "virtual" spherical body. 

In order to analyze in more details the motion of the axially symmetric 
body, we consider the following three Cartesian coordinate systems: 
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1. A principal coordinate system B attached to the axially symmetric body 
and centered at its center of mass. 

2. A coordinate system S attached to the "virtual" spherical body and cen- 
tered at its center of mass. 

3. An inertially fixed coordinate system TV. 

Assumptions 1. Without loosing generality we assume that the coordinate 
system B has its third axis parallel to the axis of symmetry of the body, and 
that the coordinate systems B and S have superimposed axes at the initial 
time (t = 0), i.e. that Rgg(0) is an identity matrix. 

We now take advantage of the reduction theorem by subdividing the 
procedure to reach the final goal of solving the motion of the coordinate 
system B with respect to the inertial frame into two subproblems: first we 
solve for the motion of the coordinate system B with respect to 5, and then 
we look for a solution of the motion of the coordinate system S with respect 
to TV. The advantage of this procedure stays in the fact that the problem of 
solving the motion for an axially symmetric body is reduced to the simpler 
problem of solving the motion for a body having spherical ellipsoid of inertia. 

Because of the definitions of cu h and uj_ e (see Eq 22 and Eq. 23) and of 
the Assumptions 1, the coordinate system B is moving with respect to the 
coordinate system S by rotating about its third axis with the angular rate 



Indeed, the vector uj h is the angular velocity of the coordinate system S with 
respect to the coordinate system TV. 

In particular, the solution for the kinematic description of the motion of 
B with respect to S can be stated as follows (see also Eq. 4 and 6) 



= |wj = A{iU_ h • g). 



(26) 



R B s(t) = RBs(0)exJ f 



cos [f(t)} sin [/(*)] 
-sin [/(*)] cos [f(t)\ 











1 



(27) 
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with 



/(*) 



(28) 



In developing Eq.27 it has been taken into account that i?ss(0), without loss 
of generality, is an identity matrix because of the Assumptions 1, and 





o > , n ( s u S b) (t) 

-w e {t) 



LU e (t) 

-w e (t) 




(29) 



In order to fully solve for the motion of the axially symmetric body it now 
remains to find the solution for the motion of the coordinate system S with 
respect to the coordinate system N. 

In particular, the dynamics of the 'virtual" spherical body with respect 
to the inertial frame is governed by Eq. 25. By expressing all of the vectors of 
that equation in scalar components along the coordinate system S, it yields 



1 



1 



R S b m, 



(30) 



I I 

where s uJh is the column matrix obtained by projecting the angular velocity 
vector Lo h along the coordinate system S, defined as 



(31) 



and B m is the column matrix obtained by projecting the torque vector m 
along the coordinate system B. 

The initial conditions for Eq. 30 are (see Eq. 24 and Assumptions 1) 



W(0) 




R S B(0) B h o 



Po 

<lo 

h r 



(32) 



where po, qo and ro are the initial conditions for the angular rate of the rigid 
body with revolution ellipsoid of inertia. 
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Furthermore, from Eq. 26 it yields 

w e (t) = Ar(t). (33) 

In particular, for all of the cases when the external acting torque depends at 
most on time, from Equations 30 it results 

u e (Jt)=A (jro + jJ Q m e {0d^j, (34) 

where m e = m ■ e. 

The kinematic equation in terms of motion of the coordinate system 5* 
with respect to N is 

Rns = Rns O ( s u h ) . (35) 

Finally, once obtained the solution for the Equations 34, 27, 30 and 35, the 
complete solution for the dynamic description of the motion of the axially 
symmetric rigid body is given at each time instant t by 

B Lu BN {t) = B cj BS (t) + R BS {t) s ui SN (t), (36) 

while the solution of the correspondent kinematic problem is 

RBN(t)=R B s{t)R S N(t). (37) 

Indeed, Eq. 36 is immediately obtained from Eq. 37 by taking the time 
derivative of both sides of Eq. 37, by considering Equations 4 and 5, and 
by exploiting the property 

ft (Rbs s ^>sn) = Rbs ft ( S usn) Rsb- (38) 



5 Exact analytic solutions for the motion of a rigid body with a 
revolution ellipsoid of inertia 

By leveraging the developments of previous sections, this section introduces 
the exact analytic solutions for the dynamics and kinematics of a rigid body 
having ellipsoid of inertia of revolution. Four cases are considered in the 
following four subsections. 
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5.1 Case with a constant external torque parallel to the symmetry axis 

The following theorem introduces a previously unknown exact analytic solu- 
tion for the complete dynamic and kinematic problems of an axially symmet- 
ric rigid body with a constant external torque staying at any time parallel 
to the symmetry axis, and arbitrary initial orientation and angular velocity. 

Theorem 2 Given poj?Oi r o and U3 real numbers, the solution of the dy- 
namic problem of determining, at any time t, the absolute angular velocity 
of a rigid body having two equal principal moments of inertia about the prin- 
cipal body axes (a\) and (0-2), initial angular velocity components p(0) = po , 
q(0) = go; an d r(0) = r$ along the three principal axes, and subjected to a 
constant torque U3, normalized by the value of the equal principal moment of 
inertia (I) and directed along the unequal inertia principal axis (03), is the 
following 



p cos [f{t)\ + q sin [/(f)] 
-Po sin [/(f)] + q cos [/(f)] 
r + t 3 U ^ 



(39) 



while the solution of the correspondent kinematic problem of determining, at 
any time t, the orientation of the body with respect to the inertial frame, i.e. 
of the coordinate system B with respect to N , is given by 



R B N{t) = R BS {t)RsN{t) (40) 

where RsN{t) is obtained from Eq. 16 of Corollary 2, by substituting the 
subindex B by S , the Rbk matrix by the identity matrix, the values of the 
initial conditions by p' = po,q'o = qo, r' = -f ro and the value of the torque 
magnitude by U — U3. Finally, Rbs (t) is given by Eq. 27 with 



f(t) = A^fr t+^t^ 



(41) 
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Proof In this case we suppose to have the following external torque acting 
on the axially symmetric rigid body 



m(t) = I }, (42) 




with m3 a scalar constant. 

From Equations 27 and 42, it results that the external torque has the 
same scalar components as seen from both the coordinate system B and 
the coordinate system S, which is attached to the "virtual" spherical body; 
this results from the fact that the torque is directed along the rotation axis 
between the two coordinate systems. Therefore, it yields 

s m(t) = R SB (t) B m(t) = B m(t). (43) 

By taking into account the Equations 42 and 43, the solution of the 
dynamic problem for the "virtual" spherical body, i.e. of Eq. 30 with initial 
conditions 32, is immediately given by 









Po 


s u SN (t) = < 




H 






r(t) t 







(44) 



r + U 3 t 



where Ua — ma/ 1. 

The solution of the kinematic problem for the "virtual" spherical body 
correspondent to the kinematic solution of Eq.44, stated in Eq. 35 , i.e. the 
determination of the rotation matrix RsN(t), is given by the Corollary 2 (see 
Section 3 and the statement of Theorem 2). 

Now, in order to reach the goal of fully solving the motion of the axially 
symmetric rigid body with respect to the inertial frame, it just remains to 
find the solution for the motion of the coordinate system B with respect to 
the coordinate system S. 
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In particular, b ojsb is immediately given by Eq. 29 with (from Equa- 
tions 34 and 44) 

w e (t) = A r + U 3 tj . (45) 

Furthermore, the motion of the coordinate system B with respect to S is 
governed by Eq. 27 with 

/(*) = £u e (OdZ = A(c+ I fr t+ 1 ^ t 2 ) , (46) 

where C — because of Assumptions 1. 

Finally, the complete solution of the dynamic problem is given by Eq. 39, 
which follows from the Equations 36, 44, and 45, while the solution of the 
kinematic problem is given by Eq. 40. 

Theorem 2 results therefore proven. 



5.2 Cases with a constant external torque perpendicular to the symmetry 
axis 

This section introduces the exact analytic solution for the complete dynamic 
and kinematic problems of an axially symmetric rigid body with an external 
torque which is constant in magnitude and stays at any time perpendicular 
to the symmetry axis. Two cases are presented, in two different theorems. 

In particular, the following theorem introduces a previously unknown 
exact analytic solution for an axially symmetric rigid body with a constant 
external torque staying at any time perpendicular to the symmetry axis and 
rotating at a constant rate about the symmetry axis, with arbitrary initial 
angular velocity and body orientation. 

Theorem 3 Given Po,qo,ro and U real numbers, the solution of the dy- 
namic problem of determining, at any time t, the absolute angular velocity of 
a rigid body having two equal principal moments of inertia (I) about the prin- 
cipal body axes (o~\) and (02) and principal moment of inertia I3 about the 
third axis (cr 3 ), initial angular velocity components p(0) = po , q(0) = qo, and 
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r(0) = r along the three principal axes, and subjected to a constant torque 
perpendicular to the third axis, rotating about the third axis at the angular 
rate ( 13 J 1 ) ro, and having magnitude U normalized by the value of the equal 
principal moment of inertia (I), is the following 



(pv + Ut) cos [/(*)]+ g sin [/(*)] 
B u BN (t)= { -(p + Ut)sm[f(t)] + q cos[f(t)} 

To 



> • 



(47) 



while the solution of the correspondent kinematic problem of determining, at 
any time t, the orientation of the body with respect to the inertial frame, i.e. 
of the coordinate system B with respect to N , is given by 



(48) 



where RsN(t) is obtained from Eq. 16 of Corollary 2, by substituting the 
subindex B by S , the time-independent matrix Rbk by 



Rsk = 



1 
1 

-10 



the values of the initial conditions by p' = po,q' = qo,r' 
RBs(t) is given by Eq. 27 with 



/(*) 



I 



rot. 



(49) 

If ro • Finally, 
(50) 



Proof In case of external torque perpendicular to the symmetry axis of the 
body, the dynamics and kinematics of the motion of the coordinate system 
B with respect to the coordinate system S are independent of the acting 
torque. Indeed, Eq. 34, by taking into account that m e = 0, yields 



= A 



h 



I-h 



ro- 



(51) 



Then b u>sb is immediately obtained by inserting Eq. 51 into Eq. 29. Fur- 
thermore, Rssit) is given by Eq. 27 with 



/(*) =We* 



(52) 
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In order to reach the goal of fully solving the motion of the axially sym- 
metric rigid body with respect to the incrtial frame, it just remains to be 
solved the problem of finding the motion of the "virtual" spherical body with 
respect to the inertial frame, i.e. the motion of the coordinate system S with 
respect to the coordinate system N. 

Because of the assumption made in the statement of Theorem 3 that the 
external torque is rotating about the axis of symmetry at a constant angular 
rate equal to — ui ei it yields 



m cos(w e i) 
'm(t) = { —m sm(Lu e t) 




> - 



(53) 



with m a scalar constant. 

From Equations 27 and 53, it results that the external torque has the 
following expression in components along the S coordinate system 



'm(t) = R SB {t) B m{t) 




(54) 









Po + Ut j 


s u SN {t) = < 


?(*) 


H 













By taking into account the Equations 53 and 54, the solution of the dy- 
namic problem for the "virtual" spherical body, i.e. of Eq. 30, is immediately 
given by 



(55) 



where U = m/I. 

The solution of the kinematic problem for the "virtual" spherical body 
correspondent to the kinematic solution of Eq. 55 is given by Corollary 2. 

Finally, the complete solution of the dynamic problem is given by Eq. 47, 
as it follows from the Equations 36, 51, and 55, while the solution of the 
kinematic problem is given by Eq. 48. 

Theorem 3 results therefore proven. 
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The following theorem introduces the exact analytic solution for the com- 
plete dynamic and kinematic problems of an axially symmetric rigid body 
with a constant external torque fixed with the axially symmetric body and 
perpendicular to its axis of symmetry, and with the initial absolute angular 
velocity also perpendicular to the symmetry axis and arbitrary initial body 
orientation. 

Theorem 4 Given po,qo an d U real numbers, the solution of the dynamic 
problem of determining, at any time t, the absolute angular velocity of a rigid 
body having two equal principal moments of inertia about the principal body 
axes (o\) and (o- 2 ), initial angular velocity components p(0) =Po, q{ty=qo, 
andr(0)= along the three principal axes, and subjected to a constant torque 
fixed with the axially symmetric body and perpendicular to its symmetry axis 
( without loss of generality we can assume that the torque is directed as the 
first axis of the coordinate system B ), is the following 



3 UJ BN {t) 




(56) 



while the solution of the correspondent kinematic problem of determining, at 
any time t, the orientation of the body with respect to the inertial frame, i.e. 
of the coordinate system B with respect to N , is given by 

R BN {t) = R SN (t), (57) 

where RsN(t) is obtained from Eq. 16 of Corollary 2, by substituting the 
subindex B by S , the time-independent matrix Rbk by the matrix Rsk de- 
fined in Eq. J±9, and the values of the initial conditions by p' ~ Po,q' Q = 
qo,r' = 0. 

Proof In this case lo c = 0, i.e. the coordinate system B is at rest with respect 
to the coordinate system S. In other words, in this case the axial symmetric 
body moves under the effect of the acting constant torque analogously to a 
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spherical symmetric body. Therefore, the theorem statement follows directly 
from Eq. 3 and from Corollary 2. 



5.3 Case without external torques 

A solution to this classical case (Euler-Poinsot's spontaneous motion) is well 
known (see for instance Leimanis 1965). However, an original and simpler 
solution, as regards its derivation, is presented here below in terms of the 
rotation matrix. This solution is obtained by applying Hestenes' reduction 
theorem (Hestenes 1999). Hestenes 1999 reports a similar solution but he 
uses the not widely familiar notion of Spinor of rotation in order to describe 
the body orientation, while here below the solution is reported in terms of 
the rotation matrix. 

The solution developed here for both the dynamics and kinematics are 
advantageous in that they are obtained with less algebraic steps than the 
classical solution, which relies on the geometrical construction of the body- 
fixed cone and the "space-fixed cone" . 

Theorem 5 Given po, qq, ro real numbers, the solution of the dynamic prob- 
lem of determining, at any time t, the absolute angular velocity of a rigid 
body having two equal principal moments of inertia about the principal body 
axes (o~\) and (02), initial angular velocity components p(0) =Po, q(0) = qo, 
and r(0) = ro along the three principal axes, and not subjected to external 
torques is the following 



p cos (uj £ t) + q sin (w e t) 
i LUBN{t) = ^ — p sin (ui e t) + q cos (cu e t) 
ro 



(58) 



with 0J e = r (Eq. 34). Furthemore, the solution of the correspondent 
kinematic problem of determining, at any time t, the orientation of the body 
with respect to the inertial frame, i.e. of the coordinate system B with respect 
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to N, is given by 



RBN(t) — RbsW RsN(t) 



(59) 



where RsN{t) is obtained from Eq. 6, withujsN = wsat(O) as given by Eq. 32. 
Finally, Rssit) * s given by Eq. 27 with 



Proof In this case, Eq. 25 states the conservation of the angular momentum 



Therefore, the dynamic solution is straightforwardly given by following the 
development of Section 4, with the kinematic problem solved by considering 
Eq. 6. An alternative to the use of Eq. 6 is the use of the Euler's equation 
(see for instance Goldstein 1980) 



RsnW = {cos(0(t))X+ [1 - cos(0(t))] aa T - sin(0(t)) Q{a)} R SN {0), 



where .!?(•) is the skew-symmetric matrix function appearing also in Eq. 5, 
X indicates the three by three identity matrix, a indicates the unit vector of 
the Euler's axis projected in either the S or the N coordinate system, and 
4> indicate the Euler's angle of rotation, which is satisfying, in this case, the 
following simple equation 



In writing Eq. 62 it has been taken into account that, in accordance with 
Assumptions 1, Rsn(0) — Rbn(0)- Finally, the axis of rotation a can be 
immediately obtained by considering that it coincides with the direction of 
the angular velocity s ujsn(0), therefore yielding 




(60) 



h(t) = h . 



(61) 



(62) 



<Kt) = \ s cosN(0)\t. 



(63) 



s ojsn(0) 
\ s ojsn(0)\' 



(64) 
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Interestingly, the overall spontaneous motion of an axially symmetric rigid 
body having three degrees of freedom is seen as intrinsically equivalent to the 
one degree of freedom motion of the "virtual" spherical body with respect 
to the inertial frame (motion of S with respect to N). This motion captures 
the precession of the inertial symmetry axis of the physical body about the 
constant and inertially fixed absolute angular momentum vector. On the 
other hand, the motion of the coordinate system B with respect to S captures 
the relative spinning of the axially symmetric body about its axis. 

6 Conclusions 

For the first time, to the knowledge of the author, the exact analytic solutions, 
for both the kinematic and the dynamic problems, have been introduced in 
this paper for the rotational motion of a rigid body having revolution ellipsoid 
of inertia, in the following three cases 

1. External torque constant in magnitude and parallel to the symmetry axis, 
with arbitrary initial angular velocity and body orientation. 

2. External torque constant in magnitude, perpendicular to the symmetry 
axis and rotating at a constant rate about the symmetry axis, with arbi- 
trary initial angular velocity and body orientation. 

3. External torque and initial angular velocity perpendicular to the symme- 
try axis, with the torque constant in magnitude and fixed with the body 
and arbitrary initial body orientation. 

The kinematic solutions are presented in terms of the rotation matrix. 

In particular, the results of this paper are obtained by building upon the 
recently found exact analytic solution for the motion of a rigid body with a 
spherical ellipsoid of inertia, and by decomposing the motion of an axially 
symmetric rigid body into the combination of the motion of a "virtual" spher- 
ical body with respect to the inertial frame and that of the axially symmetric 
body with respect to this "virtual" body. 
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All of the analytical solutions introduced in this paper have been suc- 
cessfully checked against numerical solutions obtained through differential 
equations propagation, for sample cases. 

Longuski 1991 and Tsiotras and Longuski 1991 give an exact analytic 
solution for the dynamic equation of an axially symmetric body which is 
applicable in a more general case (torque fixed with the body and having 
a generic direction) than the cases introduced in the present paper. Never- 
theless the exact analytic solutions introduced in the present paper, when 
they apply, are critically advantageous because they are complete for both 
the dynamics and the kinematics. 

In conclusion, the newly introduced exact analytic solutions are of high 
theoretical interest, as they add new elements to the set of few special cases 
for which a complete exact solution of the rotational motion of a rigid body 
exists. Additionally, they are also of interest from the applied mathematics 
point of view as they constitute significant new comparison cases for the 
validation and error analysis of approximate algorithms. In particular, the 
analytic solutions presented in this paper are valid for any length of time and 
rotation amplitude. 
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